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Abstract— Based on the epicyclic approximation, we have simulated the motion of the 
young open star clusters IC 4665 and CoUinder 359. The separation between the cluster 
centers is shown to have been minimal 7 Myr ago, 36 pc. We have established a close 
evolutionary connection between IC 4665 and the Scorpius-Centaurus association — the 
separation between the centers of these structures was ~ 200 pc 15 Myr ago. In addition, 
the center of IC 4665 at this time was near two well-known regions of coronal gas: the Local 
Bubble and the North Polar Spur. The star HIP 86768 is shown to be one of the candidates 
for a binary (in the past) with the pulsar PSR B1929+10. At the model radial velocity of 
the pulsar K = 2 ± 50 km s~^, a close encounter of this pair occurs in the vicinity of IC 
4665 at a time of —1.1 Myr. At the same time, using currently available data for the pulsar 
B1929+10 at its model radial velocity Vr = 200 ±50 km s~^, we show that the hypothesis of 
Hoogerwerf et al. (2001) about the breakup of the COph — B1929+10 binary in the vicinity 
of Upper Scorpius (US) about 0.9 Myr ago is more plausible. 

INTRODUCTION 

At present, fairly reliable estimates of such parameters as the distance and proper motion 
components are available for a number of nearby neutron stars, but as yet there is no infor- 
mation about their radial velocities. In the problem of determining the possible birthplace 
of a neutron star, the radial velocity is a free parameter. 

Based on the Hipparcos Catalogue and radioastrometric data, Hoogerwerf et al. (2001) 
analyzed 56 nearby high- velocity stars (r < 700 pc), the so-called runaway stars, to determine 
their possible birthplace inside nearby young stellar groups — open star clusters or associ- 
ations. One of their conclusions was that the star C Oph and the pulsar PSR B1929+10 
(which, according to the J2000 coordinate system, was designated as PSR J1932+1059) 
originated jointly in a supernova explosion in a binary about 1 Myr ago in the part of the 
Scorpius-Centaurus association designated as US (Upper Scorpius). 

However, subsequent VLBI measurements of the parallax for PSR B1929+10 (Brisken 
et al. 2002; Chatterjee et al. 2004) showed that the probable distance to the pulsar differs 
significantly from that used by Hoogerwerf et al. (2001). In the opinion of Chatterjee et al. 
(2004), the new data basically reject the hypothesis of Hoogerwerf et al. (2001) about the 
birthplace of PSR B1929+10 at the boundary of the Scorpius-Centaurus association and its 
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connection with the star ( Oph. In this paper, we test Hoogerwerf 's hypothesis using new 
observational data. In addition, estimating the probability of other hypotheses is also of 
considerable interest. 

One of such hypotheses was proposed by Walter (2001), who suggested considering a 
different neutron star as a binary (in the past) with ( Oph, namely, RX J 185635-3754. 

For our part, we also suggest considering the hypothesis that the birthplace of PSR 
B1929-I-10 could be the open cluster IC 4665 or Cr 359. The idea of this hypothesis is 
based on the simulations of the PSR B1929-I-10 trajectory over the celestial sphere (Fig. 2 
from Hoogerwerf et al. (2001)) at various radial velocities (+200 or —200 km s^^), which 
showed that it passed in the immediate vicinity of the open clusters IC 4665 and Cr 359. In 
this paper, we determine the radial velocity of the pulsar at which its encounter with these 
clusters is closest. 

In the immediate solar neighborhood with a radius of 200-300 pc, we know such regions 
of interstellar tenuous hot gas with a temperature T 10^ K as the Local Bubble and the 
North Polar Spur. Since the Local Bubble is a compact region that is almost free from 
absorbing matter, it was first discovered by Fitzgerald (1968) through an analysis of the 
interstellar reddening of stars. The absorption maps of the distribution of neutral gas in Nal 
hues constructed by Sfeir et al. (1999) show an inclination to the Galactic plane as for the 
Gould Belt. An overview of the formation scenarios for the North Polar Spur can be found 
in Miller et al. (2007). Data on all of the known bubbles in the solar neighborhood were 
collected in Heiles (1998). 

The physical processes that produce bubbles, in outline, are as follows. Multiple super- 
nova explosions occur in young open clusters. This gives rise to stellar winds and shock 
waves that sweep out the gas to the periphery of the local region in question, where it 
produces blobs in the form of shells or extended walls. At the shell boundaries, the gas 
density increases significantly, the gas cools down, and molecular clouds are formed. During 
supernova explosions, the gas inside the shell can heat up and radiate in the X-ray range. In 
the opinion of Berghofer and Breitschwerdt (2002), the most reahstic theory for the origin 
of the Local Bubble is the hypothesis about multiple, but not simultaneous explosions of 
about 20 supernovae during the past 10-20 Myr. Frish (1995), Maiz-Apcllaniz (2001), and 
Breitschwerdt and de Avillez (2006) substantiate the viewpoint that supernova explosions in 
the Scorpius-Centaurus association, which is part of the Gould Belt, are the most suitable 
objects responsible for the formation of both the Local Bubble and the North Polar Spur. 

The open clusters IC 4665 and Cr 359 have a high probability of belonging to the Gould 
Belt both in their age and spatial position (Piskunov et al. 2006) and in kinematics (Bobylev 
2006). Since there are very few young open star clusters in the first Galactic quadrant that 
are members of the Gould Belt and that have rehable kinematic data (particularly the radial 
velocities), studying the trajectories of the open clusters IC 4665 and Cr 359 and members 
of the Scorpius-Centaurus association is of great importance for understanding the evolution 
of the Gould Belt. 

The goal of this paper is to establish the evolutionary connections of the young open 
clusters IC 4665 and Cr 359 with the structure of the Gould Belt, to establish the possible 
influence of the clusters on the bubble development, and to determine the most probable 
birthplace of the pulsar PSR B1929-I-10 by varying its radial velocity. 

The problems posed here are solved by constructing the trajectories of the objects using 
the epicyclic approximation. 
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THE METHOD 



In this paper, we use a rectangular Galactic coordinate system with the axes directed away 
from the observer toward the Galactic center {I = 0°,b = 0°, the X axis), along the Galactic 
rotation (/ = 90°, 6 = 0°, the Y axis), and toward the North Galactic Pole (6 = 90°, the Z 
axis) . The corresponding space velocity components of the objects U, V, W are also directed 
along the X,Y,Z axes (Kuhkovskii 1985). 

The method of the epicyclic approximation (Lindblad 1927, 1959) allows the orbits of 
stars to be constructed in a coordinate system rotating around the Galactic center in a 
circular orbit. We assume that the origin of the coordinate system coincides with the local 
standard of rest and that the stars move along epicycles in the direction opposite to the 
Galactic rotation. We use the method in a form associated with a rectangular coordinate 
system (Fuchs et al. 2006): 

X{t) = X{0) + ^sm{Kt) + ^(1 - cos(Kt)), (1) 

Y{t) = Y{Q) + 2a(x{Q) + - ^ViO) smiKt) + ^t/(0)(l - cos(Kt)), 

Z(t) = ^^^^ sm(ut) + Z(0) cos(ut), 

where t is the time in Myr (pc/Myr=0.978 km s"^), which is measured into the past; k = 
\/— 4f2o-B is the epicyclic frequency; A and B are the Oort constants; Qq is the angular 
velocity of Galactic rotation of the local standard of rest, Qq — A — B; and v = ^/inGpo is 
the frequency of the vertical oscillations, where G is the gravitational constant and po is the 
star density in the solar neighborhood. The space velocities of the objects are calculated for 
any necessary time using the formulas 

U(t) = U(0) cos(Kt) ^V^(O) sin(Kt), (2) 

V{t) = ^^C/(0) sm{Kt) + 1/(0) cos(Ki), 

K 

W{t) = W{0)cos{iyt) - Z{0)usm{ut). 

The parameters X{0),Y{0), Z{0) and ^(0), V{0),W{0) in Eqs. (1) and (2) denote the current 
positions and velocities of the objects. The velocities U, V, W are given relative to the local 
standard of rest with {U, V, W)lsr = (10.00, 5.25, 7.17) ± (0.36, 0.62, 0.38) km s"^ (Dehnen 
and Binney 1998). 

Following Fuchs et al. (2006), we adopted po = 0.1 MqPC~^, which gives u = 0.074 km 
s~^ pc~^. We also took the Oort constants A = 13.7 ± 0.6 km s"^ kpc^^ and B = —12.9 ± 
0.4 km s~^ kpc~^ that we found previously (Bobylev 2004) by analyzing the independent 
determinations of these parameters by various authors; in that case, k — 0.037 km s~^ pc~^. 

THE DATA 

The necessary input data on the objects under consideration, such as their equatorial coor- 
dinates, proper motion components, radial velocities, and parallaxes, are given in Table 1. 



Table 2 gives the spatial coordinates -'^(O), 1^(0), Z{0) and space velocities U{0), V{0), W{0) 
of the stars and clusters under consideration. 

IC 4665 and Cr 359. For the open cluster IC 4665, we use the coordinates of the center, 
proper motion components, and an estimate of the photometric distance r = 352 ± 70 pc 
from the CRVOCA Catalog (Kharchenko et al. 2007) and the radial velocity from Manzi 
et al. (2007), where it was determined using 39 most probable cluster members. According 
to Manzi et al. (2007), the age of IC 4665 is 27.7l3;5 Myr- This estimate is interesting in 
that it depends little on the quality of isochrones, because it was obtained by comparing 
the lithium abundances in cluster stars with such well-known young clusters as NGC 2547 
and IC 2391. Other authors give the following "isochronic" cluster age estimates: 36 Myr 
(Mermilliod 1981), 30-100 Myr (Prosser 1993), and 43 Myr (Piskunov ct al. 2006). 

For the open cluster Cr 359, we use the coordinates of the center, proper motion com- 
ponents, and radial velocity from the CRVOCA Catalog (Kharchenko et al. 2007) and the 
distance r = 450 ± 200 pc from Lodieu et al. (2006), who conducted a critical review of the 
estimates obtained by various authors. According to Lodieu ct al. (2006), the age of Cr 359 
is 60 ± 20 Myr. Other authors estimate this cluster as a younger one, 32 Myr (Piskunov et 
al. 2006) and ^ 30 Myr (Wielen 1971; Abt and Cardona 1983). 

According to Kharchenko et al. (2007), the apparent radii of IC 4665 and Cr 359 are 
l.°0 and l.°l, respectively. 

de Wit et al. (2006) estimated the mass of IC 4665 to be 300-350 Mq. Cr 359 probably 
has a similar mass. 

PSR B1929+10. The radio pulsar PSR B1929+10 is an isolated neutron star with 
an age of 3 Myr (estimated from the ratio P/2P and belongs to the population of nearby 
neutron stars that are closely related to the Gould Belt (Popov et al. 2003; Motch et al. 
2006). 

For PSR B1929-I-10, we use the data from Chatterjee et al. (2004) that were obtained 
from VLBI measurements. The most important refined parameter is the new parallax of the 
pulsar, TT = 2.77 ± 0.07 mas. 

Note that Hoogcrwcrf ct al. (2001) adopted a model value of vr = 4 ± 2 mas. The 
preliminary distance estimates for the pulsar that made it closer to the Sun served as a basis 
for this choice. 

US, UCL, and LCC. The initial positions of the centers (X(0), y(0), Z(0)) = 
(134,-20,52) pc and heliocentric velocities {U{0),V{0),W{0)) = (-6.7,-16.0,-8.0) ± 
(5.9,3.5,2.7) km 8"^ for US, (119,-67,31) pc and (-6.8,-19.3,-5.7) ± (4.6,4.7,2.5) km 
s-i for UCL, (62, -100, 10) pc and (-8.2, -18.6, -6.4) ± (5.1, 7.3, 2.6) km s^^ for LCC were 
taken from Fernandez et al. (2006). These authors obtained the positions of the centers (with 
errors of 1.2 pc) using data from de Zeeuw et al. (1999); the velocities were determined by 
Sartori et al. (2003) based on the space velocities of a large number of stars (more than 120 
for each of the groups). Fernandez ct al. (2006) used these data to analyze the kinematics 
of the members of the Scorpius. Ccntaurus association. However, the orbits of the objects 
were calculated using a model that, apart from the axisymmetric Galactic potential, also 
included the potential from a spiral density wave and the potential of the bar at the Galactic 
center (Fernandez et al. 2008). Comparing the results obtained by various methods is of 
considerable interest. 
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RESULTS 



The Neighborhoods of IC 4665, Cr 359, and PSR B1929+10 

The encounters of the centers. Based on Eqs. (1) and (2), we obtained several solutions 
at various radial velocities of the pulsar PSR B1929+10, Vr (PSR). The results are the 
following: 

(1) at Vr (PSR) = +45 km s~^, the pulsar encounters with the cluster IC 4665 at the 
distance = V^X^ + AY^ + AZ"^ = 52 pc at t = -0.85 Myr; 

(2) at Vr (PSR) = —60 km s~^, the minimum distance from the pulsar to the center of 
Cr 359 is = 38 pc at i = -1.0 Myr; 

(3) at Vr (PSR) = km s~^, we obtain a joint encounter with the two clusters. In 
this case, the minimum distances from the pulsar to the centers of Cr 359 and IC 4665 are 
Ar = 69 pc and A^ = 60 pc, respectively, at t = — 1.0 Myr. 

As we see, the pulsar trajectory can pass through the coronas of the two clusters under 
consideration. As was noted by de Wit et al. (2006), the tidal radius of IC 4665 is 1°. 
At an assumed distance to the cluster of Or = 352 pc, it is ~ 6 pc. However, the stars with 
a common proper motion occupy an area of ^ 100 square degrees on the celestial sphere 
(Lodieu et al. 2006) and, hence, the coronal radius of IC 4665 is ^ 30 pc. As can be seen 
from Table 2, the error in the spatial position of the center of IC 4665 is about ±50 pc; its 
spatial displacement in 2 Myr is insignificant. 

In this respect, Cr 359 is of lesser interest, since the error in its distance is 44%. Therefore, 
the radius of its corona lies within the range 20-60 pc, but the spatial localization accuracy 
of the cluster is low, about ^ 150 pc (Table 2). 

Next, we take the next natural step or, more specifically, we search for a suitable rapidly 
flying star, a possible member of a binary. For this purpose, we use the list of such stars 
from Hoogcrwcrf et al. (2001). We can see from Fig. 2 in Hoogcrwcrf et al. (2001) that 
the trajectories of three stars passed in the immediate vicinity of IC 4665 and Cr 359 on the 
celestial sphere with a radius of ^ 15° 2 Myr ago: HIP 66524, HIP 86768, and HIP 91599. 
To calculate the space velocities, we invoke the current radial velocities of these stars from 
the Pulkovo Catalog of Radial Velocities (Gontcharov 2006). 

The construction of the trajectories for the three rapidly flying stars shows that HIP 
86768 is most interesting among them. 

Indeed, HIP 66524 moves away fairly rapidly toward the Galactic center in X coordinate, 
without having a close encounter with the objects of interest to us. 

For HIP 91599 (star No. 20 from the fist of Hoogerwerf et al. 2001), the situation is 
slightly different. At the radial velocity of the pulsar PSR B1929+10 Vr (PSR) =+250 km 
s~^, it encounters with HIP 91599 at the distance A,. = 47 pc at t = —0.5 Myr. The modulus 
of the velocity difference between the pulsar and HIP 91599 at the time of their encounter 
is Av = VAU^ + AV^ + AW^ = 316 km s-\ 

At the radial velocity of the pulsar Vr (PSR)=+2 km s~^, it encounters with HIP 86768 
at the distance A^ = 19 pc at i = —1.1 Myr. In this case, the modulus of the velocity 
difference between the pulsar and HIP 86768 at the time of their encounter Ay = 206 km 
s^^. Taking into account the fairly large errors in the distances and velocities of HIP 86768, 
we can assume that a closer encounter is also possible. Therefore, simulating the encounters 
by taking into account the errors in the data is of great interest. 



Monte Carlo Simulations of Encounters. (1) Initially, we directly repeated the 
numerical experiment of Hoogerwerf et al. (2001) on the encounter between the star ( Oph 
and PSR B1929+10. We computed 3 miUion orbits by taking into account the random errors 
in the input data, which were normally distributed with the 3a region. For the pulsar, we 
used old data: VriPSR) = 200 ± 50 km s'^ tt = 4 ± 2 mas, /i„ cos5 = 99 ± 12 (6 x 2) mas 
yr~^, and /is = 39 ± 8 (4x2) mas yr~^ (following Hoogerwerf et al. (2001), we doubled the 
initial errors in the proper motions, as shown in parentheses). The results are: out of the 3 
miUion orbits, 32505 encounters occur at distances < 10 pc (1.1%); in 4410 of the 32505 
cases, the star and the pulsar were no farther than 10 pc from the US center about 1 Myr 
ago. The derived parameters are in excellent agreement with the results of Hoogerwerf et 
al. (2001). 

(2) The same experiment on the encounter between ( Oph and PSR B 1929+10 with the 
currently available data for the pulsar: K (PSR) = 200 ± 50 km s-\ tt = 2.77 ± 0.7(0.07 x 
10) mas, fiaCosS = 94.09±3.3(0.11 x 30) mas yr-^ and i^s = 42.99±4.8(0.16 x 30) mas yr'^ 
The initial values of the random errors (as shown in parentheses) were increased in order to 
have conditions comparable to those in experiment (1). The results are: out of the 3 million 
orbits, 74115 encounters occur at distances < 10 pc (2.5%); in 5611 of the 74115 cases, 
the star and the pulsar were no father than 10 pc from the US center about 1 Myr ago. 
The domains of admissible values of Vr, tt, HaCOsS (/i^), and fj.s at which 5611 encounters 
occur are presented in Fig. 1 for the star C Oph and in Fig. 2 for PSR B1929+10. Figure 3 
shows the expected distribution F3D of minimum distance A^ calculated from the formula 
of Hoogerwerf et al. (2001), 
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(3) 



for the adopted a = 2.5 pc. 

(3) The experiment on the encounter between HIP 86768 and PSR B1929+10 for Vr 
(PSR) = 2 ± 50 km s^^ and the same values of tt = 2.77 ± 0.7(0.07 x 10) mas, /i„ cos 5 = 
94.09 ± 3.3(0.11 X 30) mas yr'^ and fis = 42.99 ± 4.8(0.16 x 30) mas yr^^ as those in 
experiment (2). The results are: out of the 3 million orbits, 22708 encounters occur at 
distances A^ < 10 pc (0.8%); in 6932 of the 22708 cases, the star and the pulsar were no 
father than 80 pc from the center of IC 4665 about 1 Myr ago. The domains of admissible 
values of V,., tt, fi^cosd (/U*), and fig at which 6932 encounters occur are presented in Fig. 4 
for HIP 86768 and in Fig. 5 for PSR B1929+10. Figure 6 shows the expected distribution 
of minimum distance A^ calculated from Eq. (3) for the adopted a — 6 pc. 

(4) The experiment of the simulation of encounters between HIP 91599 and PSR 
B1929+10 at K (PSR) = 250 ± 50 km s^^ shows that out of the 3 million orbits, 32961 
encounters occur in the range 15 pc< A^. < 30 pc (0 encounters at A^. < 15 pc) about 0.5 
Myr ago. The stars HIP 86768 and HIP 91599 do not encounter at distances A^ < 10 pc, 
suggesting that the triple system HIP 86768— HIP 91599— B1929+10 is unlikely. 

(5) Testing the hypothesis of Walter (2001) — the experiment on the encounter between 
C Oph and the isolated neutron star (NS) RX J185635-3754 at K (NS)= -50 ± 50 km 
s^^. The results are: out of the 3 million orbits, 2144 encounters occur at distances A^. < 
10 pc (0.07%); there are no encounters with the US center at A^ < 10 pc. The derived 
characteristics are in good agreement with the results of Hoogerwerf et al. (2001). 
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IC 4665 and the Scorpius-Centaurus Association 

We traced the trajectories of IC 4665 and Cr 359 in a time interval comparable to their 
lifetime, up to t = —30 Myr. For the mean values of the input parameters described in 
Section 2.1, we found the minimum separation between the cluster centers to be = 36 pc 
Bit t — —7 Myr. Allowance for the errors in the cluster distances and velocities showed that 
the encounter time lies within the interval 0-12 Myr. 

Figure 7 shows the positions of IC 4665, Cr 359, and members of the Scorpius-Centaurus 
association and their trajectories in the past 30 Myr; since the age of US does not exceed 5 
Myr, the corresponding part of the trajectory is marked by the dotted line. 

Figure 8 shows the calculated positions of IC 4665, the members of the Scorpius- 
Centaurus association UCL and LCC, and two bubbles, the Local Bubble (LB) and the 
North Polar Spur (LI) (Breitschwerdt and de Avillez 2006) at t = —15 Myr. The difference 
in the positions of the UCL and LCC centers at this time found from comparison with the 
data of Ortega et al. (2002) or Fernandez et al. (2006) is ±50 pc along the X and Y 
axes and very small, ±3 pc, along the Z axis. We think these discrepancies attributable 
to differences in the models to be insignificant for the purposes of this paper. Note the 
approach of Fuchs ct al. (2006), who simulated not the trajectories of the group centers, but 
the individual trajectories of very massive stars, potential supernovae, using Eq. (1). Some 
of them turned out to be near bubble LI at t = —15 Myr (see Fig. 8). For IC 4665, apart 
from the coordinates of the center. Fig. 8 shows a cloud of 300 points (distributed with the 
3cr region) calculated by taking into account the errors in the input data for the cluster. As 
can be seen from the figure, the contribution from the errors in the distance is dominant. 

DISCUSSION 

The Connection of IC 4665 and Cr 359 with the Scorpius-Centaurus 
Complex 

As we see from Fig. 7, the edge of the Gould Belt located in the first Galactic quadrant 
continuously approached the Galactic plane in a time interval of ^ 30 Myr in the past. This 
is in good agreement with the model for the evolution of the Gould Belt suggested by Olano 
(2001) and is consistent with the model calculations of the kinematic evolution of the Gould 
Belt performed by Perrot and Grenier (2003) using data on molecular clouds for the scenario 
with a "Galactic-plane crossing" (Fig. 7 from Perrot and Grenier (2003)), whereby the age 
of the Gould Belt is 51.8 ± 1.0 Myr. According to these two models, the edge of the Gould 
Belt located in the first Galactic quadrant reaches the Galactic plane at ~ — 10 Myr. Our 
results show a similar motion. Thus, for example, the Z coordinates of the centers of UCL 
and LCC (which lie in the plane of the Gould Belt) at t = — 15 Myr fixed in Fig. 7 are close 
to pc, while the Z coordinates of the centers of IC 4665 and Cr 359 (rising above the plane 
of the Gould Belt) are close to 70 pc. 

Most of the authors believe that the Scorpius-Centaurus association is part of the Gould 
Belt (de Zeeuw et al. 1999; Bobylev 2006). However, Fernandez et al. (2008) advocate the 
point of view that the association could evolve independently of the Gould Belt (outside the 
scope of the supernova explosion hypothesis) and could be formed from the parent cloud 
that was compressed by collisions with a spiral density wave. 
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On the other hand, whatever the origin of the association, the history of star formation 
in it agrees satisfactorily with the model of successive star formation suggested by Blaauw 
(1964; 1991) and developed by Preibish and Zinnecker (1999), as applied to US. The age 
estimates for the members of the Scorpius-Ccntaurus association (de Geus et al. 1989) lie 
within the ranges 5-6 Myr for US, 14-15 Myr for UCL, and 11-2 Myr for LCC. Both the 
present-day "isochronic" estimates (Mamajek et al. 2002; Sartori et al. 2003), 8-10 Myr 
for US and 16-20 Myr for UCL and LCC, and our kinematic age estimate for the entire 
association, 21 Myr (Bobylev and Bajkova 2007), agree well with these estimates. 

Our results show that the separation between IC 4665 and the Scorpius-Centaurus as- 
sociation was much smaller in the past. Thus, for example, whereas the current separation 
between their centers is 302 pc (UCL-IC 4665), it was 215 pc 15 Myr ago and 120 pc 30 
Myr ago. The minimum separation was ~ 80 pc at —66 Myr, but our model approximation 
on such a long time scale is already unrehable. All of these facts suggest that they were 
formed from the same parent hydrogen cloud; the cluster IC 4665 has always been located 
on its periphery. 

According to the model of successive star formation, the shocks from supernova explosions 
compress the nearby (within 100 pc) molecular clouds, triggering star formation in them. 
In our view, the role of IC 4665 in the evolution of the Scorpius-Centaurus association could 
lie in the fact that supernova explosions in this cluster could trigger star formation processes 
in the association itself. 

At t = —15 Myr, IC 4665 was near two well-known regions of coronal gas: the Local 
Bubble and the North Polar Spur. 

Based on an analysis of HI data, Heiles (1998) obtained the following parameters for 
the North Polar Spur (which is known as the brightest component of radio loop I): the 
heliocentric distance to the center is 120 pc, the coordinates of the center are / = 320° and 
b — +5°, and the radius is 118 pc. According to the model of Willingale et al. (2003), it can 
be represented as an expanding spherical super bubble with the following parameters: the 
distance is 210 pc, the direction of the center is I — 352°, b — +15°, and the bubble radius is 
140 pc. X-ray observations of the North Polar Spur showed that the peak of emission with 
an energy of 3/4 keV is observed in the direction with / = 26. °8 and b = +22. °0 (Snowden 
et al. 1997; Wilhngale et al. 2003; Miller et al. 2007), i.e., at the bubble boundary. The 
current coordinates of IC 4665 are I — 30. °6 and b — -|-17.°1. Taking into account the results 
shown in Figs. 7 and 8, we conclude that IC 4665 has always been located at a distance 
of 150-200 pc from the bubble boundary both in the past and at present. This means that 
supernova explosions in IC 4665 over the past several Myr could influence the formation 
of the North Polar Spur. This influence lies in the fact that supernova explosions in IC 
4665 could produce a stellar counter-wind (the main explosions in the Scorpius-Centaurus 
association), which additionally compresses the bubble walls. 

The Probable Birthplace of the Pulsar PSR B 1929+ 10 

We considered the hypothesis of Hoogerwerf et al. (2001) that the star ( Oph and the pulsar 
B1929-I-10 could be components of a binary in the neighborhood (within less than 10 pc) 
of US about 1 Myr ago using currently available data for the pulsar. Our simulations of 
encounters showed that, contrary to the view of Chatterjee et al. (2004), improved data for 
the pulsar only strengthen this hypothesis. Therefore, it seems most plausible. This requires 
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that the pulsar radial velocity fall within the range Vr = 200 ± 50 km s~^. 

Our hypothesis of the HIP 86768 — B1929+10 binary, which could exist about 1 Myr 
ago in a wide neighborhood (within less than 80 pc) of IC4665, also seems plausible. The 
probability of the existence of such a binary is lower than that in the hypothesis of Hoogerwerf 
et al. (2001). A distinctly different pulsar radial velocity, 2 ± 50 km s~^, is required 

for our hypothesis to hold. 

An examination of the hypothesis by Walter (2001) about the ( Oph — RXJ185635-3754 
binary showed its probability to be an order of magnitude lower than that of our hypothesis. 



CONCLUSIONS 

Based on the cpicyclic approximation, we simulated the motion of the young open star 
clusters IC 4665 and Collinder 359. The separation between the cluster centers was found 
to have been minimal 7 Myr ago, 36 pc. This suggests that IC 4665 and CoUinder 359 were 
formed from the same parent hydrogen cloud. 

Wc showed a close evolutionary connection between IC 4665 and the Scorpius-Ccntaurus 
association. Thus, for example, the separation between the centers of these structures was 
~ 200 pc about 15 Myr ago. At that time, IC 4665 was near two well-known regions 
of coronal gas: the Local Bubble and the North Polar Spur. This means that supernova 
explosions in IC 4665 over the past 15 Myr could affect the development of bubbles and, in 
particular, the North Polar Spur. 

Analysis of the parameters for the encounter of the pulsar PSR B1929+10 with various 
stars in the vicinity of IC 4665 and Cr 359 suggests that the star HIP 86768 is a suitable 
candidate for a binary (in the past) with the pulsar. Monte Carlo simulations of the en- 
counters of this pair with the pulsar radial velocity K = 2 ± 50 km s'^ showed that out of 
the 3 million orbits, 22708 encounters occur at distances < 10 pc (0.8%); in 6932 of the 
22708 cases, the star and the pulsar were no farther than 80 pc from the center of IC 4665 
about 1 Myr ago. Wc showed that using currently available data for the pulsar increases 
the probability of the hypothesis by Hoogerwerf et ai. (2001) about the breakup of the C 
Oph — B1929-I-10 binary in the immediate vicinity of US about 1 Myr ago. Monte Carlo 
simulations of the encounters of this pair with the pulsar radial velocity Vr — 200 ± 50 km 
s~^ showed that out of the 3 million orbits, 74115 encounters occur at distances A,. < 10 pc 
(2.5%); in 5611 of the 74115 cases, the star and the pulsar were no farther than 10 pc from 
the US center about 1 Myr ago. 
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Table 1: Data on the objects 



Object 


Ct( J2000.0) , 
r 

0(J2000.0) 


/ia COS 6, 
— 1 

mas yr 


IJ'S, 

—1 

mas yr 


mas 


Vr, 
1 „— 1 

km s 


B1929+10 


19''32"*13^9497 
10 59 32 .4203 


94.09 ±0.11 


42.99 ±0.16 


2.77 ±0.07 




IC 4665 


5 43 


-0.57 ±0.30 


-7.40 ±0.36 


2.84 ±0.56 


-16.0 ± 1.1 


Cr 359 




0.22 ±0.28 


-8.90 ±0.26 


2.22 ±0.99 


-4.6 ±0.2 


HIP 66524 


13''38™09^0095 
-50°20'59".702 


-60.75 ± 0.49 


-13.34 ±0.67 


2.28 ±0.75 


23 ± 1 


HIP 81377 
(C Oph) 


16"37™09^5378 
-10°34'01".524 


13.07 ±0.85 


25.44 ±0.72 


7.12 ±0.71 


-9.9 ±5.5 


HIP 86768 


17M3"'47^0205 
-7°04'46".588 


-7.12 ±0.70 


-10.39 ±0.51 


2.34 ±0.80 


19.0 ±4.3 


HIP 91599 


18^0^48^0517 
-8°43'07".688 


-9.64 ± 1.13 


-22.64 ±0.79 


3.61 ± 1.16 


29 ±5 


RXJ185635-3754 


18'*56'"35^56 
-37°54'37".0 


326.7 ±0.8 


-59.1 ±0.7 


16.5 ±2.3 





Table 2: Positions and heliocentric velocities of the objects 





X(0), 


Y{0), 


Z(0), 


^(0), 


V{0), 


W{0), 




pc 


pc 


pc 


km s ^ 


km s ^ 


km s~^ 


B1929+10 


244 ±6 


265 ±7 


-25 ± 1 


-109 ±34 


91 ±37 


-107 ±6 


IC 4665 


290 ± 57 


171 ± 34 


104 ± 21 


-6 ±2 


-17 ±3 


-9 ±2 


Cr 359 


381 ± 169 


217 ±96 


99 ±44 


6±8 


-16± 14 


-10 ±8 


HIP 66524 


279 ± 92 


-326 ± 107 


90 ±30 


-83 ±67 


-102 ±57 


1 ±3 


HIP 81377 


128 ± 13 


14 ± 1 


56 ±6 


-11 ±5 


18 ±2 


-1±2 


HIP 86768 


397 ± 136 


134 ± 46 


86 ±29 


25 ±7 


-18± 17 


6±2 


HIP 91599 


253 ±81 


113 ±36 


-8 ±3 


39 ±7 


-18± 13 


-3± 1 


RXJ18B63B-3764 


58 ±8 


-1.4 ±0.2 


-18 ±3 


-75 ± 48 


20 ± 1 


-75 ± 15 



Note. The velocities of the pulsar B 1929+ 10 were calculated for its radial velocity = 2 ± 50 km 
s~^; the velocities of the pulsar RX J185635.3754were calculated for Vr = —50 ± 50 km s~^. 
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Fig. 1. Domains of admissible values at which 5611 encounters occur at distances < 10 pc 
between the star C Oph and the pulsar PSR B1929+10, during which they fall into the 
neighborhood of US, < 10 pc for ( Oph. 



400 




LOO 200 300 4DD 




1 1 


J , ■ 1 


1 1 





100 200 300 40) 



Fig. 2. Domains of admissible values at which 5611 encounters occur at distances < 10 pc 
between the star C Oph and the pulsar PSR B1929+10, during which they fall into the 
neighborhood of US, < 10 pc for PSR B1929+10. 



14 




Fig. 3. (a) Expected distribution of minimum distance < 10 pc for 5611 encounters of 
( Oph with PSR B1929+10 and (b) histogram of encounter times. 
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Fig. 4. Domain of admissible values at which 6932 encounters occur at distances A,. < 10 pc 
between HIP 86768 and PSR B1929+10, during which they fall into the neighborhood of IC 
4665, < 80 pc for HIP 86768. 
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Fig. 5. Domains of admissible values at which 6932 encounters occur at distances < 10 pc 
between HIP 86768 and PSR B1929+10, during which they fall into the neighborhood of IC 
4665, < 80 pc for PSR B1929+10. 
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Fig. 6. (a) Expected distribution of minimum distance A,. < 10 pc for 6932 encounters of 
HIP 86768 with PSR B1929+10 and (b) histogram of encounter times. 
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Fig. 7. Spatial coordinates of IC 4665, Cr 359, members of the Scorpius-Centaurus associa- 
tion, and their trajectories in the past 30 Myr. 
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Fig. 8. Positions of IC 4665, UCL, LCC, the Local Bubble (LB), and the North Polar Spur 
(LI) (Breitschwerdt and de Avillez 2006) 15 Myr ago.. 
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